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The dynamics of a soliton propagating in a single-mode optical fiber with gain, loss, 
and Raman coupling to thermal phonons is analyzed. Using both soliton perturbation 
theory and exact numerical techniques, we predict that intrinsic thermal quantum 
noise from the phonon reservoirs is a larger source of jitter and other perturbations 
CN ! than the gain-related Gordon-Haus noise, for short pulses (< lps), assuming typical 

fiber parameters. The size of the Raman timing jitter is evaluated for both bright and 
dark (topological) solitons, and is larger for bright solitons. Because Raman thermal 
quantum noise is a nonlinear, multiplicative noise source, these effects are stronger 
S^h. for the more intense pulses needed to propagate as solitons in the short-pulse regime. 

Thus Raman noise may place additional limitations on fiber-optical communications 



Abstract 



and networking using ultrafast (subpicosecond) pulses. 
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1. Introduction 



In this paper, we analyze in some detail the effects of Raman noise on solitons. In particular 
we derive approximate analytic expressions and provide further detail for the precise nu- 
merical results published earlieiS. The motivation for this study is essentially that coupling 
to phonons is one property of a solid medium that definitely does not obey the nonlinear 
Schrodinger equation. The presence of Raman interactions plays a major role in perturbing 
the fundamental soliton behaviour of the nonlinear Schrodinger equation in optical fibers. 
This perturbation is in addition to the more straightforward gain/loss effects that produce 
the well-known Gordon-Haus effect!. 

The complete derivation of the quantum theory for optical fibers is given in an ear- 
lier paper!, denoted (QNI). That paper presented a detailed derivation of the quantum 
Hamiltonian, and included quantum noise effects due to nonlinearities, gain, loss, Raman 
reservoirs and Brillouin scattering. Phase-space techniques allowed the quantum Heisenberg 
equations of motion to be mapped onto stochastic partial differential equations. The result 
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was a generalized nonlinear Schrodinger equation, which can be solved numerically or with 
perturbative analytical techniques. 

The starting point for this paper is the phase-space equation for the case of a single 
polarization mode, obtained using a truncated Wigner representation, which is accurate in 
the limit of large photon number. We use both soliton perturbation theory and numerical 
integration of the phase-space equation to calculate effects on soliton propagation of all 
known quantum noise sources, with good agreement between the two methods. 

Our main result is that the Raman noise due to thermal phonon reservoirs is strongly 
dependent on both temperature and pulse intensity. At room temperature, this means that 
Raman jitter and phase noise become steadily more important as the pulse intensity is in- 
creased, which occurs when a shorter soliton pulse is required for a given fiber dispersion. 
Using typical fiber parameters, we estimate that Raman-induced jitter is more important 
than the well-known Gordon-Haus jitter for pulses shorter than about one picosecond. Al- 
though we do not analyze this in detail here, we note that similar perturbations may occur 
during the collision of short pulses in a frequency- multiplexed environment. 

2. Raman- Schrodinger Model 

We begin with the Raman-modified stochastic nonlinear Schrodinger equation [Eq. (6.3) of 
(QNI)], obtained using the Wigner representation, for simplicity: 
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Here = tyyvt /n is a dimensionless photon field amplitude, while r = (t — x/v)/t and 
C = x/xo, where t is a typical pulse duration used for scaling purposes and x = tl/\k"\ is 
a characteristic dispersion length. The group velocity v and the dispersion relation k" are 
calculated at the carrier frequency loq. 

Apart from a cut-off dependent vacuum noise, the photon flux is J = \<f)\ 2 n/t , where 
n = \k"\Ac/(n 2 hujQt ) = v 2 t /xxo is the typical number of photons in a soliton pulse of 
width to, again for scaling purposes. In this definition, the fibre is assumed to have a modal 
cross-sectional area A and a change in refractive index per unit intensity of n 2 . The positive 
sign in front of the second derivative term applies for anomalous dispersion (k" < 0), and 
the negative sign applies for normal dispersion (k" > 0). The functions g and h are gain/loss 
and Raman scattering response functions respectively, while T and T R are stochastic terms, 
discussed below. 

Similar, but more accurate, equations occur with the positive-P representation, although 
in this case, the phase-space dimension is doubled. In order to simplify the calculations 
further, we assume that gain and loss in the fiber are broadband relative to the soliton 
bandwidth, and balance exactly. This requires that the amplifier sections in the fiber are 
sufficiently close together (of the order of the soliton scaling length or less) so that the soliton 
can propagate without distortion!. 
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For the analytic calculations, we also assume that the Raman nonlinear response function 
is instantaneous on the timescale of the soliton width. This is equivalent to assuming that 
the phonon modes are heavily damped, and means that the Raman coupling leads to only 
incoherent scattering of the propagating radiation. While this approximation neglects the 
well-known self-frequency shiftpnD, we find that the self-frequency shift by itself is not a major 
cause of jitter for the distance scales we consider here. This assumption can be improved at 
the expense of more complicated analytic calculations. However, the full equations are used 
in the numerical simulations, which agree quite well with our analytic predictions. 

The Raman-modified equation then reduces to 

0(r,C) + r c (r,C) , (2.2) 

where the bracketed term represents the usual nonlinear Schrodinger (NLS) equation in 
normalized, propagative form. The combined noise sources have been grouped together as 

r c (r,C)=r(r,C)+^(r,C)0(r,C). (2.3) 

A. Initial conditions and quantum evolution 

Equation (|2.2|) is a complex-number equation that can accurately represent quantum oper- 
ator evolution through the inclusion of various noise sources. In the absence of any noise 
sources, this equation reduces to the classical nonlinear Schrodinger equation. This deter- 
ministic limit corresponds to taking n — > oo. As well as the noise sources explicitly appearing 
in Eq. ( |2.2|) , there must be noise in the initial conditions to properly represent a quantum 
state in the Wigner representation. Regardless of the initial quantum state chosen, there 
must be at least minimal level of initial fluctuations in to satisfy Heisenberg's uncertainty 
principle. We choose to begin with a multimode coherent state, which contains this minimal 
level of initial quantum noise and which is an accurate model of mode-locked laser output. 
This is also the simplest model for the output of mode-locked lasers, and we note that, 
in general, there could be extra technical noise. For coherent inputs, the Wigner vacuum 
fluctuations are Gaussian, and are correlated as 

(A0(r, O)A0*(r', 0)) = ±_5(r - r'). (2.4) 

Physical quantities can be calculated from this phase-space simulation by averaging prod- 
ucts of (f> and 0* over many stochastic trajectories. In this Wigner representation, these 
stochastic averages correspond to the ensemble averages of symmetrically-ordered products 
of quantum operators, such as those representing homodyne measurements and other mea- 
surements of phase. 

B. Wigner noise 

Both fiber loss and the presence of a gain medium each contribute quantum noise to the equa- 
tions in the symmetrically-ordered Wigner representation. The complex gain/absorption 
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noise enters the nonlinear Schodinger (NLS) equation through an additive stochastic term 
T, whose correlations are: 

(r(n, or {of, CO) = {aG ^ n aA) 5(c - Om + no, (2.5) 

where T(Q, () is the Fourier transform of the noise source: 

1 r°° 

T{n,C) = -7= rfrr(r,C)exp(^r). (2.6) 



. 2tt 

The dimensionless intensity gain and loss are given by a G and a A , respectively. 

Similarly, the real Raman noise, which appears as a multiplicative stochastic variable 
T R , has correlations 



(T R (Q, C)T R (Q', CO) = =5(C - C')W + nO 

n 



n th (Q) + - 



(2.7) 



where the thermal Bose distribution is given by n t h(Q) = [exp (fr|f2| //c#Tt ) — 1]~ and where 
a R (Q) is the Raman gain, whose profile is given in Fig. 1 of (QNI). Thus the Raman noise 
is strongly temperature dependent, but it also contains a spontaneous component which 
provides vacuum fluctuations even at T = 0. 

As the n dependence of all the noise correlations show, the classical limit of these quantum 
calculations is the deterministic nonlinear Schrodinger equation. The problem of jitter in 
soliton communications is an example of how intrinsic quantum features can have a direct 
macroscopic consequence, even in a way that impinges on current developments of applied 
technology. There are, of course, classical contributions to jitter, such as noise arising from 
technical sources. However, it is the jitter contributions from essentially quantum processes, 
namely spontaneous emission in fibre amplifiers, that is the current limiting factor in soliton 
based communications systems. Other jitter calculations rely on a classical formulation with 
an empirical addition of amplifier noise, and important predictions of the Gordon-Haus effect 
have been obtained. Nevertheless, this quantum treatment presented here of all known noise 
sources is necessary to determine the limiting effects of other intrinsic noise sources, which 
become important for shorter pulses and longer dispersion lengths. 

In the absence of the noise sources, the phase-space equations have stationary solutions 
in the form of bright (+) or dark (— ) solitons. Solitons are solitary waves in which the 
effects of dispersion are balanced by nonlinear effects, to produce a stationary pulse that is 
robust in the presence of perturbations. We note here that, in reality, the Raman response 
function is noninstantaneous, which causes a redshift in the soliton frequency. This soliton 
self- frequency shift is a deterministic effect, and so can be neglected in the treatment of noise 
effects, to a first approximation. The accuracy of this approximation will be evident in the 
subsequent comparison of analytic with numerical results. The numerical results all include 
the complete nonlinear response function, rather than the approximate instantaneous form 
given above. 

Excessive self-frequency shift may cause problems when finite bandwidth elements are 
used. However, it has been shown! that bandwidth- limited gain can in fact cancel the effect 
of the Raman redshift, by pulling the soliton back towards the centre of the spectral band. In 
the simulations we show in this paper, the total redshift is estimated to be Af ~ 0.02THz, 
which is small compared to the total width of the gain spectrum in typical fiber laser 
amplifiers (Av ~ 3THz)i. 
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3. Perturbation Theory 

We now proceed to derive the approximate analytic expressions for the effects of noise on 
soliton jitter, using soliton perturbation theoryErfl, for both bright and dark solitons. 



A. Bright solitons 



The stationary soliton of Eq. fl2.2| ) for anomalous dispersion is: 
^bright( r > = Asech[Ar - q(()} expfiVr + 10(C)], 



(3.1) 



where dq/dQ = VA and 86 /d( = (A 2 — Y21L2, with amplitude A and velocity V. Fol- 
lowing the method presented by Haus et at^, we treat the effects of the noise terms as 
perturbations around a soliton solution whose parameters vary slowly with (: 



0(r ) C) = 0(r,C) + A0(r,C), 

where the unperturbed soliton solution is given by: 
0(r, C) = A(C)sech[A(C)r - q(()) exp\iV{C)r + i9{$] 



for a bright soliton. Substituting Eq. ( |3.2| ) into Eq. 
equation (first order in A0(r, ()) 



(3.2) 
(3.3) 

gives the following linearized 
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A0(r, C) 



±^ + *20*(r,O0(r,C) 



A0(r, C) + i(f>(r, C) 2 A0*(r, C) + T(r, C) , (3.4) 



where the linearized noise source T(r, Q is defined as: 

r(r l c) = r(r,c)+<r B (r l c)?(r,c) • 



(3.5) 



Now we wish to determine the evolution of the soliton parameters as a function of 
propagation distance (. To do this, we expand the perturbation in terms of the soliton 
parameters plus a continuum term: 



i 1 

= £/ Pi AP i + A0 c (T,C). 

i 

where Pi G {V, q, A, 6}. The projection functions for each parameter are 
1 



(3.6) 



A 



— r tanh(v4r — q) 



f q = tanh(Ar — 
f v = ir$, 
fe = i~4>- 



(3.7) 
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Since the linearised equation Eq. (|3.4j) is not self-adjoint, these eigenfunctions are not or- 
thogonal. In order to select out the evolution of particular parameters, we therefore choose 
an alternative set of functions: 

h = 

fq = T~4>, 

fv = itanh(Ar — q)(p, 

fo = iTta,nh(AT-q)0. (3.8) 



These are the eigenfunctions of the adjoint equation Eq. (|3.4| ), and obey the orthogonality 
condition 

* {II' 1 ' fr-fr/} kr (3.9) 

Substituting the Taylor expansion [Eq. (|3.6|) 1 into the linearized equation [Eq. ( |3.4|) 1 and 
using the functions fp i to project out particular parameters shows that the growth of fluc- 
tuations in position Aq is governed by 

^Ag(C) = AAV(C) + r fl (C) 

^A^(C) = MC), (3.10) 

where we have taken the unperturbed velocity to be zero: V = 0. The stochastic terms are 
defined as 

r^cc) =»{jT drV(0F(r,C)}. (3-11) 

Here we have assumed that the perturbations in the continuum <p c are orthogonal to the 
fp i . This depends on such perturbations dispersing sufficiently rapidly away from the region 
around the soliton. In fact, any nonsoliton perturbation will disperse and would also move 
away from the soliton, since the group velocity for any linear perturbations will be different 
to the propagation velocity of the soliton. 

We wish to find the growth of fluctuations in position q((). Because the position de- 
pends on the soliton frequency V, the contributions arising from both Y q and Ty must be 
considered. Firstly, 

r g (C) = ^{/°° dTATsech(AT - q) exp(-iVr - i0)r(r,C)} 

/oo 
dTATsech(Ar - q)dt{exp(-iVr - iO)T}, (3.12) 
-oo 

and 

IV(C) =^y°° dTA(-i)sech(AT - q) tanh(Ar - q) exp(-tVr - i6)T{r, j 

/oo _ , 

drAsech(AT - q) tanh(Ar - q) Asech(Ar - q)T R + ^{exp(-iVr - i6)T} . (3.13) 
-oo L - 1 
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From this we can calculate the growth of the fluctuations in velocity: 
10 

| jH A0(r, ()fv*dr} + £ dCTy(C') 



AV(C)=AV(0)+ [ C d('T v (C) 
Jo 



(3.14) 



Using the noise correlations calculated above, the correlations in the velocity fluctuations 
can now be calculated: 



(av(c)av(0> = (av(o)au*(o)>+ / c ^V^nMC'OrMO) 

JO JO 



i4 
6n 



+ 



Q! -A | 2A 2 X(t ) 

3n 



n 



c c < c, 



(3.15) 



where the overlap integral I (to) is defined as 

POO /'OO _ 

X(to)= / drdr'tanh(r)sech 2 (r)tanh(r')sech 2 (r / )J : '(r/A-r7A). 



(3.16) 



— oo «/ — oo 



Here .F(t) is the inverse Fourier transform of the fluorescence jF(f2) = -[/^(fl) + |]a ii (fi). 

The correlations in position fluctuations correspond to the jitter in arrival times, be- 
cause we have chosen a propagative reference frame. The jitter therefore feeds off position 
fluctuations as well as noise entering through the velocity: 



(A g (C)Ag*(CO) = (Ag(0)Ag*(0)> 

+ / c f C d("dC" \a 2 (av((")av*(C")) + (r g (C")r*(C ,// )) 

Jo Jo 

Thus the timing jitter is 
<[Ar(C)] 2 > = (Ag(C)Ag*(C)) 



(3.17) 
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a G A 3 2A i l(t ) 



9n 



+ 



3n 



c 3 



(3.18) 



which contains cubic terms due to the gain and Raman couplings, and also slower growing 
terms due to the initial vacuum fluctuations and amplifier noise. 

We note that an alternative method that exploits conserved quantities in the NLS equa- 
tion is often use for deriving the timing jitter. The linearised approach that we have 
presented has the advantage that derivatives of products of stochastic variables do not ap- 
pear. With such derivatives, the normal rules of calculus do not apply. Rather, the rules of 
Ito stochastic calculus must be observed, leading to extra drift terms. 

B. Dark solitons 



Fibers in the normal dispersion regime can support dark soliton solutions, so called since 
they correspond to a dip in the background intensity^: 
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ark( r ' = <W 1 - -4 2 sech [(f) Ar - q(()] exp[i0(C)] exp[ia(C, r)], 
a(C, r) = arcsm i V , (3.19) 



1 - A 2 sech 2 [0 o ^-g(C)] 

where d9/d( = 0q, dq/d( = Ay/1 — A 2 <b\ and O is the amplitude of the background field. 
The size of the intensity dip at the center of the soliton is given by A, with the intensity 
going to zero in a black tanh(r) soliton, for which A = 1. Dark solitons are classed as 
topological solitons, because they connect two background pulses of different phase. The 
total phase difference between the boundaries is ip = 2arcsin(A). 

The nonvanishing boundary conditions of the dark pulse complicate the perturbation 
calculation of jitter variance. To ensure that all relevant integrals take on finite values, we 
impose periodic boundary conditions at r = ±77, which are taken to infinity at the end of 
the calculation. These boundary conditions require a soliton solution of the form 

{ib ib ib I 

cos - + i sin - tanh[0 o r sin - - g(C)j > , (3.20) 

with a wavenumber offset k = ~ arctan tan ^ tanh(0 o ^ sin ^) . The perturbation theory 
now proceeds in a similar fashion to the bright soliton case, except that we can greatly 
simplify the calculation if the unperturbed solution is taken to be a black soliton, i.e. ib = 71". 
The projection functions for the soliton parameters Pi G {ib, q, 0o, 9} are 

fe = — 0o tanh(0 O ' r — q) exp(i9 — inr) 
f<f> = «[tanh(0 o r — q) + o rsech 2 (0o r — q)] exp(i9 — ikt) 
f g = — z0osech 2 (0 O T — q) exp(i9 — ikt) 

U = 0o[/^i0o^tanh(0 o r - q) - -) exp(i9 - ikt), (3.21) 
where f3\ = l/[20o r / tanh(0 O T/)]. For the required adjoint functions, we choose: 
f g = — ^-^sech 2 (0 O T — q) exp(i9 — ikt) 

U = n 1 sech 2 (0 o ^ - q) exp(i9 - ikt), (3.22) 
— Pi — 1 

where 7g = 4/(3 J^ n dtsechH) and 7^ = (ft - 1)/ J^ n dt(ftt tanhtsech 2 t-sech 2 t/2). The 
orthogonality condition is now 

/ dTf P J P ; = S itj . (3.23) 

J -n+q/4>o J 

Once again, the adjoint functions can be used in the linearized equation [Eq. (p.4|)j to 
determine how the fluctuations in position evolve: 

|A g(C ) = ^f A^(C) + r 9 (C) 

^A^(C) = IV(C), (3.24) 



s 



where /?2 = 4>o tanh(0 O ' r «) — l/ r z- Here we see how the fluctuations in phase produce fluctu- 
ations in position. The stochastic term in the equation for ip evaluates to 



IV (C) = & 



drU*(C)r(r,0 



dr 



... , sech 2 (0oT - q) ${exp(-i6 + ikt)Y} - O tanh(0 o r - q)? , (3.25) 

-n+q/4> Pi — 1 L J 

from which the correlations of the phase fluctuations can be calculated: 

K rC 



(A^(C)A^(CO) = (A^(0)A^*(0)> + 

2 r 2a u Tv > 



7v> 



+ 



o jo 

GU, 2 



«'<iV(C")r;(0> 



+ 



27^,% (t c 



3n 7g (pi - l) 2 0o n(A - 1) 



3n 7g o 

where the overlap integral X n (to) is now defined as 



c c < c, 



drdr' tanh(r)sech 2 (r) tanh(r / )sech 2 (r')jF(r/0 o — t'/^o 
The leading order terms for the fluctuations in position are thus 



^3.26) 



(3.27) 



(Ag(C)Ag*(C)) 



4 Jo Jo 
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12n7g 



^ 2 + 



«'(A?KC")A^(C)> 

« G 0o/3l 7y> 
18rr 7g (pi 



6n 



c 3 . 



(3.28) 



By taking the limit 77 — > 00, we find the leading order terms in the jitter growth for a black 
soliton to be 



<[Ar(C)] 2 > 



12n 



°c 2 + 



18n 6n 



(3.29) 



where the overlap integral Z(£o) is as defined in Eq. (|3.16|) . As in the anomalous dispersion 
regime, the vacuum fluctuations contribute to quadratic growth in the jitter variance, and 
gain and Raman fluctuations contribute to cubic growth. However, the size of the jitter 
is smaller than that in the bright soliton case, for the same propagation distance (. The 
contribution from the vacuum and gain terms is one half and the contribution from the 
Raman term is one quarter of that in Eq. ( |3.18 ), giving dark solitons some advantage over 
their bright cousins. 



4. Scaling properties 

In summary, there are three different sources of noise in the soliton, all of which must be 
taken into account for small pulse widths. These noise sources contribute to fluctuations in 
the velocity parameter, which lead to quadratic or cubic growth in the timing-jitter variance 
for single-pulse propagation. The noise sources also produce other effects, such as those 
effected through soliton interactions, but we will not consider these here. 

Each of the noise sources has different characteristic scaling properties, which are sum- 
marized as follows: 
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A. Vacuum Fluctuations 



The vacuum fluctuations cause diffusion in position which is important for small propagation 
distances. There are position fluctuations even at the initial position, since the shot noise 
in the arrival time of individual coherent-state photons gives an initial fluctuation effect. 
After propagation has started, this initial position fluctuation is increased by the additional 
variance in the soliton velocity, due essentially to randomness in the frequency domain. 
For bright solitons the resulting soliton timing variance is given byM 

<[MC)] 2 >/ = ^= + ^C 2 (bright). (4.1) 

For purposes of comparison, note that N = 2n is the mean photon number for a sech(r) 
soliton. Numerical calculations confirm that for tanh(r) dark solitons, the variance was 
about one half the bright-soliton value, as predicted by the analysis in outlined in Sec. 



^^■ = im + ^ (dark) - (4 - 2) 

This shot-noise effect, which occurs without amplification, is simply due to the initial 
quantum-mechanical uncertainty in the position and momentum of the soliton. Because 
of the Heisenberg uncertainty principle, the soliton momentum and position cannot be spec- 
ified exactly. This effect dominates the Gordon-Haus effect over propagation distances less 
than a gain length. However, for short pulses, this distance can still correspond to many 
dispersion lengths - thus generating large position jitter. We note that there are also initial 
fluctuations in the background continuum, which may feed into the soliton parameters as the 
soliton propagates. This lessor effect is included in the numerical calculation; a comparison 
of the numerical results with the analytic results confirms that the initial fluctuations in the 
soliton parameters account for almost all of the shot-noise contribution to the jitter. 



B. Gordon-Haus noise 

As is well known, the noise due to gain and loss in the fiber produces the Gordon-Haus 
effect, which is currently considered the major limiting factor in any long-distance soliton- 
based communications system using relatively long (> lOps) pulses. Amplification with 
mean intensity gain a G , chosen to compensate fiber loss, produces a diffusion (or jitter) in 
position. Unless other measures are taken, for sufficiently small amplifier spacingp and at 
large distances this is given by§B0 

<[Ar(C)] 2 > G „ ~ — C 3 (bright), 

rv 

([Ar(C)] 2 )^ ^ ^=C 3 (dark), (4.3) 

in which the linearly growing terms have been neglected. 

Another effect of the amplifier noise is to introduce an extra noise term via the fluctua- 
tions in the Raman-induced soliton self-frequency shift. This term scales as the fifth power 



10 



of distance and hence will become important for long propagation distances. This com- 
bined effect of spontaneous emission noise and the Raman intrapulse scattering has been 
dealt with by others!. The full phase-space equation [Eq. (|2.1|) 1 models this accurately, 
since it includes the delayed Raman nonlinearity, and the effect would be seen in numerical 
simulations carried out over long propagation distances. 



C. Raman noise 

A lesser known effect are the fluctuations in velocity that arise from the Raman phase-noise 
term T R in Eq. ( |2.2|) . Like the Gordon-Haus effect, this Raman noise generates a cubic 
growth in jitter variance: 

<[Ar(C)] 2 > R = (bright), 

<[Ar(C)] 2 > R = ^C 3 (dark). (4.4) 



where X(to) is the integral defined in Eq. (|3.16| ) that indicates the spectral overlap between 



the pulse spectrum and the Raman fluorescence. The mean-square Raman induced timing 
jitter has a cubic growth in both cases, but the dark soliton variance is one quarter of that 
of the bright soliton. 

The magnitude of this Raman jitter can be found by evaluating X(t ) numerically, or 
else using an analytic approximation. An accurate model of the Raman gain, on which 
I(t ) depends, requires a multi-Lorentzian fit to the experimentally measured spectrumta. 
A fit with 11 Lorentzians was used in the numerical simulations, including 10 Lorentzians 
to accurately model the measured gain and fluorescence. One extra Lorentzian was used at 
low frequencies to model GAWBS (Guided Wave Acoustic Brillouin Scattering); this has a 
relatively small effect on an isolated soliton, except to causephase noise. 

For analytic work, however, a single-Lorentzian modeled can suffice for approximate 
calculations. A plot of the Raman gain profile a R (Q) for both models is given in (QNI), along 
with a table of the fitting parameters for the multi-Lorentzian model. The spectral features of 
the Raman noise correlations are determined directly from the Raman fluorescence function 
•F(O), which we plot in Fig. [I]. For the single-Lorentzian model, the fluorescence spectrum 
is approximately flat at low frequencies: 



1 
2 

^2 



n th (Q) + ]- 



2F 1 n 1 6jk B T = 



which greatly simplifies the Raman correlations. As Fig. [I] indicates, the spectral overlap 
of jF(f2) with a to = lps soliton occurs in this low frequency region. Thus the white noise 
approximation for the Raman correlations is good for solitons of this pulse width and larger. 
For smaller pulse widths, not only is the Raman contribution to the noise larger due to the 
greater overlap, but the colored nature of the correlations must be taken into account. 
In the single-Lorentzian model, X(t — ► oo) ~ -^JF(O), which gives 
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<[At(x)] 2 > fl * ' ' 2 , ° 1 V = -JL^i ( - ) (bright), 



45Aito 45n \x / 

2|fc // | 2 n 2 ^gj zr (0) 3 _ 2^(0) /^\ 3 
45^4c*g ' 45rl \x ) 



([AtWn. B ^ 1 ' "V = If) (dark). (4.6) 



At a temperature of 300if, .F(O) = 4.6 x 10~ 2 when a single Lorentzian centered at 12 THz 
with fitting parameters Fx = 0.7263, 5 X = 20 x 10 12 t and fi x = 75.4 x 10 12 t is used. 



5. Numerical results 

More precise results can be obtained by numerically integrating the original Wigner phase- 
space equation, Eq. ( |2.1| ), which includes the full time-delayed nonlinear Raman response 
function. The results for t = 500/s bright and dark solitons are shown in Figs. |^(a) 
and 0(b), respectively. The gain and photon number were chosen to be G = a G /xq = 
4.6 x 10 _5 m _1 (0.2dB/km) and n = 4 x 10 6 , with xq = 440m. These values are based on A = 
40[/im] 2 , k" = 0.57[ps] 2 /km, and n 2 = 2.6 x lO^^fm] 2 /!^ for a dispersion-shifted fibei@. 
These numerical calculations use the multiple-Lorentzian model of the Raman response 
function shown in Fig. |l|, which accurately represents the detailed experimental response 
function. 

The numerical method is based on the split-step ideai as adapted to Raman 
propagationil. Noise is treated using a central difference technique appropriate to stochastic 
equations!!, with the necessary adaptations required to treat a partial stochastic differential 
equation^. All calculations were duplicated using two different space steps, but with the 
same underlying noise sources, in order to calculate discretization error. Sampling error was 
also estimated using standard central limit theorem procedure over a large ensemble of noise 
sources. Tests on time steps and window sizes were also carried out to ensure there were no 
errors from these sources. 

The initial conditions consist of a coherent laser pulse injected into the fiber. In the 
Wigner representation, this minimum uncertainty state leads to the initial vacuum fluctu- 
ations. The numerical calculation thereby includes the full effect of these zero-point fluc- 
tuations, including the noise that appears in the background continuum and in the soliton 
parameters. We take the initial pulse shape in the anomalous dispersion regime to be a 
fundamental bright soliton, with A = 1 and V = 9 = q = 0. Such a soliton can be experi- 
mentally realized with a sufficiently intense pulse, which will reshape into a soliton or soliton 
train. The nonsoliton part of the wave will disperse and any extra solitons will move away 
at different velocities from the fundamental soliton. The numerical simulation in the normal 
dispersion regime used two black solitons of opposite phase chirp (A = ±1), so that the field 
amplitude at either boundary was the same. This phase matching ensured the stability of 
the numerical algorithm, which assumes periodic boundary conditions. 

The position jitter at a given propagation time was calculated by combining the wave- 
form with a phase-matched local-oscillator pulse that had a linear chirp in amplitude, and 
integrating the result to give the soliton position. This homodyne measurement involves 
symmetrically ordered products, so the Wigner representation will give the correct statis- 
tics. The variance in soliton position was then calculated from a sample of 1000 trajectories. 
For this small distance of propagation (~ 10km), the jitter variance due to the initial noise 
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is twice the Gordon-Haus jitter, but for larger distances, the cubic effects are expected to 
dominate. 

For ultrafast pulses, the Raman jitter dominates the Gordon-Haus jitter (by a factor of 
two in the 500/s bright soliton case) and will continue to do so even for long propagation 
distances. For short propagation distances, the Gordon-Haus effect is not exactly cubic, 
because of neglected terms in the perturbation expansion, which give a linear (as opposed to 
cubic) growth in the jitter variance. However, there are no such terms in the Raman case. 
The analytic Raman results are also shown in the figures, which show that our approximate 
formula gives a reasonable fit to the numerical data even for subpicosecond pulses. Using 
this approximate formula, the relative size of the two effects scales as 

(bright), 

(dark), (5.1) 
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These equations show why experiments to dateE9, which have used longer pulses (to > lps) 
and dispersion-shifted fiber, have not detected the Raman-noise contribution to the jitter. 
The Raman jitter exceeds the Gordon-Haus jitter for bright solitons with periods x$ < 
1.5km. Dark solitons, on the other hand, have an enhanced resistance to the Raman noise, 
which means that a shorter period is needed before the Raman jitter will become important. 

The total jitter, which corresponds to the realistic case where all three noise sources are 
active, is also shown in Fig. [| and, in the bright soliton case, is about a factor of three larger 
than the ordinary Gordon-Haus effect, over the propagation distance shown. The physical 
origin of these quantum noise sources cannot easily be suppressed. The initial vacuum- 
induced timing jitter is caused by the shot-noise variance in the soliton guiding frequency. 
The physical origin of the Raman jitter is in frequency shifts due to soliton phase modulation 
by the ever-present quantum and thermal phonon fields in the fiber medium. 

The numerical method should give accurate results far beyond the distance shown in 
Fig. |2|, provided that the transverse and propagative resolutions are made large enough. 
The equations generated in the Wigner method should remain valid up to ( ~ \/^> which 
corresponds to about 1000km. When the Wigner equations can no longer be trusted, the 
positive- P equations will still give accurate results. This paper has not analysed any multi- 
soliton effects, although the numerical method does simulate interactions between solitons. 
This just requires the initial conditions and simulation window width to be set up accord- 
ing to whether interactions are to be considered or not. We have not included third-order 
dispersion in our model, which would become important for very small pulses (t Q ~ 100/s), 
but it could easily be included into the equations for numerical simulation. 

The approximate analytic results are most limited probably by their exclusion of Raman 
intrapulse effects, such as the deterministic self- frequency shift and the amplifier jitter that 
feeds through this. Approximate calculations^ of the self-frequency shift jitter variance show 
that it grows as the fifth power of distance. With our parameters and the measured value 
of the Raman time constantly, it would become larger than the usual Gordon-Haus effect at 
x ~ 100km, or about 10 times the propagation distance shown in Fig. |2]. Using standard 
techniques^, the perturbation theory presented in this paper could be extended to include 
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the self-frequency shift contributions (from both the amplifier noise and Raman phase noise) 
to the total jitter. 

6. Conclusions 

Our major conclusion is that quantum noise effects due to the intrinsic finite-temperature 
phonon reservoirs are a dominant source of fluctuations in phase and arrival time, for sub- 
picosecond solitons. For longer solitons, Raman effects are reduced when compared to the 
Gordon-Haus jitter from the laser gain medium that is needed to compensate for losses. 
The reason for this is the smaller intensity of the pulse, and therefore the reduced Raman 
couplings that occur for longer solitons - which are less intense than shorter solitons with 
the same dispersion. The ratio can be calculated simply from the product Gx , which gives 
the gain per soliton length. A smaller x corresponds to a shorter, more intense soliton and 
hence a larger Raman noise; while a larger G corresponds to increased laser gain, with larger 
spontaneous noise. 

At a given pulse duration and fiber length, a strategy for testing this prediction would 
be to use short pulses with dispersion-shifted fiber having an increased dispersion, since this 
increases the relative size of the Raman jitter. The physical reason for this is very simple. 
Solitons have an intensity which increases with dispersion if everything else is unchanged. At 
the same time, the multiplicative phase noise found in Raman propagation is proportional 
to intensity, and hence becomes relatively large compared to the additive Gordon-Haus 
noise due to amplification. For large enough dispersion, the temperature-dependent Raman 
jitter should become readily observable at short enough distances so that amplification is 
unnecessary. This would give a completely unambiguous signature of the effect we have 
calculated. A mode-locked fiber soliton laser would be a suitable pulse source, owing to the 
very short (64/s), quiet pulses§!il that are obtainable. 
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FIGURES 



Fig. L Spectrum of the fluorescence function !F[oj) for the 11-Lorentzian model (continuous 
lines) and the single-Lorentzian model (dashed lines), for a temperature of T = 300-fT. Also shown 
is the spectrum of a to = lps soliton. 

Fig. 2. Timing jitter in to = 500/s bright (a) and dark (b) solitons due to initial quantum 
fluctuations (circles), Gordon-Haus effect (crosses) and Raman noise (plus signs). The asterisks 
give the total jitter and the continuous line gives the approximate analytic results for the Raman 
jitter. 
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